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ABSTRACT 

A procedure for applying covariance analysis to de- 
termine the most efficient estimation strategy for 
satisfying the stringent mission requirements of long 
arc orbit determination of applications satellites is 
presented. The procedure is applied to the problem 
of satisfying mission requirements with respect to 
altitude determination of GEOS-C. It is shown that 
requirements are met when twelve dominant geo- 
potential coefficients are estimated along with satel- 
lite state. This application of covariance analysis 
is general and can be applied to future applications 
satellites . Recommendations for future studies are 
also given. 
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ESTIMATION STRATEGIES FOR ORBIT DETERMINATION 
OF APPLICATIONS SATELLITES 


INTRODUCTION 

The use of satellite based observations for earth and ocean physics applications 
has resulted in unprecedented demands for orbit determination accuracy. For 
example, the GEOS-C satellite, scheduled for launch in late 1974, will be 
equipped with an altimeter capable of one meter precision. In order to make 
effective use of such altimeter data it is advisable to determine the GEOS-C alti- 
tude with one meter or better precision. Other applications such as marine 
position determination and navigation and polar motion determination also call 
for high accuracy orbit determination. See (1) for details. 

This paper presents a procedure for utilizing covariance analysis and propaga- 
tion techniques to obtain optimal estimation strategies for long arc orbit deter- 
mination. The procedure is applied to the problem of long arc orbit determina- 
tion of GEOS-C and it is shown that the resultant accuracies are consistent with 
the needs of an earth and ocean physics program. These results are significant 
in themselves with respect to the requirements of the GEOS-C mission and they 
are also significant in that they suggest what proper estimation strategies can 
accomplish in determining orbits of future applications satellites. 
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COVARIANCE ANALYSIS AND OPTIMAL ESTIMATION STRATEGIES FOR 
LONG ARC ORBIT DETERMINATION 

It is convenient to think of the errors in an orbit determination as arising from 
two distinguishable sources. First, the data which is processed in the orbit 
determination procedure is corrupted with high frequency noise. The second 
error source is due to a misrepresentation of parameters in the functional re- 
lationship between the estimated state and the data. These parameters fall into 
two categories— measurement parameters and dynamic parameters. Measure- 
ment parameters are involved in defining the geometrical relationship between 
the satellite at a given instant and an observation of the satellite obtained at that 
instant. The most significant measurement parameters are observation biases 
and station location coordinates. Dynamic parameters are used in defining the 
relationship between the state of the satellite at one instant and the state of the 
satellite at a later instant. The most commonly mentioned dynamic parameters 
are those which determine atmospheric drag and solar pressure and the coeffi- 
cients of the spherical harmonic expansion of the geopotential field. The stand- 
ard deviation of the estimate of a given component of satellite state is the root 
sum square of a contribution due to data noise, a contribution due to uncertain 
measurement parameters, and a contribution due to uncertain dynamic parame- 
ters. To obtain an intelligent estimation strategy for determining a given orbit 
it is necessary to know the major contributions to orbit estimation error. If 
for instance the dominant error source for a given orbit and a given tracking 
configuration were due to data noise, one could recommend the use of a more 
accurate data type or the addition of extra tracking stations as the effective way 
to improve the orbital estimate. If, however, the dominant error source were 
due to uncertain parameters the addition of more data would be of little use. In 
fact, the error contribution due to data noise is seldom significant in orbit de- 
termination problems. Generally for short arc orbit determinations, uncertain 
measurement parameters provide the dominant error source and for long arc 
orbit determinations, uncertain dynamic parameters are the major contributors. 
See, for instance (2) for a study of this phenomenon. 

IIow can one determine which misrepresented parameters are significant sources 
of error for a specific orbit determination problem ? The significance of a 
given parameter as an error source, is a function both of its analytical impact 
on the functional relationship between the data and the estimated state and of 
the extent to which the parameter is likely to be misrepresented. Both factors 
must be taken into account. This can be done by simulations in the following 
manner: First simulate data for a specific orbit determination problem using 
the nominal value of the parameter under investigation. Next, process the data 
in an orbit determination program taking care that the same model used to simu- 
late the data is used in the program with the exception that the given parameter 
be perturbed by what is believed to be the standard deviation of its estimate. 
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The resultant estimate is then propagated and compared at various time points 
with the assumed true orbit. The result is a time history of the effect of a one 
sigma perturbation of the parameters on the estimate of satellite state. Ob- 
viously this procedure can be very expensive. Covariance techniques provide a 
better method of obtaining the same end. The difference in approach between a 
simulation study and a covariance analysis can be described as follows: in a 
simulation, data are generated and a least squares adjustment process is ac- 
tually performed. The estimated state is then compared to an assumed true 
state and conclusions are drawn. In a covariance analysis mode, the least 
squares adjustment process is postulated rather than actually performed and 
only its associated covariance matrix is computed. Sophisticated covariance 
analysis software permits one to list at any point in an orbit the contribution to 
the error in each component of state due to each uncertain parameter and at a 
small fraction of the cost that would be incurred in obtaining the same informa- 
tion by means of simulations. The mathematical details of this procedures are 
provided in the appendix. 

Granted that covariance techniques are effective in identifying significant error 
sources in an orbit determination process , what is to be done once this informa- 
tion is obtained? The straightforward approach is to place in an estimation 
mode along with satellite state all those parameters whose uncertainties are 
significant error sources. Another procedure which is often more sparing in 
terms of computer time is the so-called "lumped parameter approach. " With 
this approach one invents an essentially fictitious set of parameters whose com- 
bined error signature in the data is similar to that which is left by the combined 
effect of the authentic misrepresented parameters. If this is possible then by 
placing the fictitious parameters in an estimation mode the errors produced by 
the authentic misrepresented parameters will be "absorbed" by the fictitious 
parameters and a good orbit determination will be the result. And if the ficti- 
tious parameters are much less in number than the authentic parameters much 
is gained in terms of computational efficiency. See (3), (4) and (5) for inter- 
esting applications of this technique. 

Without doubt the dominant error source in long arc orbit determination of an 
applications type satellite will be the misrepresentation of geopotential coeffi- 
cients. To find a small number of fictitious parameters whose combined error 
signature resembles that left by misrepresented geopotential terms is far from 
a trivial task. Consequently the authors have opted for the straightforward ap- 
proach to long arc orbit determination in which they utilize covariance techni- 
ques to identify those geopotential terms whose uncertainties are significant 
error sources in the orbit determination process and then place these terms in 
an estimation mode along with satellite state. Properly weighted a priori esti- 
mates of the estimated geopotential coefficients are assumed in the process in 
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order to avoid the normal matrix ill-conditioning which otherwise may well re- 
sult when several geopotential coefficients are estimated from tracking data as- 
sociated with a single satellite. In the next section the procedure is utilized in 
developing estimation strategies for long arc orbit determination of the GEOS-C 
satellite . 
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estimation strategies for long arc orbit 

DETERMINATION OF GEOS-C 

The results of an error analysis of GEOS-C long arc and short arc orbit deter- 
mination are reported in (2) . In that study it was determined that in the long 
arc case the only significant errors are due to uncertainties in geopotential co- 
efficients. Hence the most direct way to improve the orbital estimate is to place 
selected geopotential coefficients in an estimation mode along with satellite 
state. To determine the most efficient use of computer time we first arrange 
the geopotential terms in the order of their significance as error sources. Then 
we investigate the resultant orbit determination accuracy , say , when the first 
three terms in the list are estimated, the first six terms, etc. 

The data postulated for the covariance analysis is collected over a seven day 
span from nine laser stations. The stations were located at Antigua, Bermuda, 
Canal Zone, Cape Canaveral, Grand Turk, Goddard, Rosman, Athens, and 
Carnarvon. The first seven of these stations are located in the Caribbean re- 
gion. For several passes per day tracking geometry is excellent in this region. 
Figure 1 demonstrates this geometry for one such pass. Two such passes are 
obtained in successive orbits. In the next six orbits visibility is poor and then 
the pattern is repeated. Since the GEOS-C period is 102 min, this implies a 
strong 13 hour observability pattern. This is the explanation for the fact that 
from both simulations and covariance analysis the orbital errors are periodic 
through the entire 7 days tracking interval with period approximately 13 hours 
(2). Consequently it is adequate to view the orbital errors during any 13 hour 
interval in order to understand the behavior of the errors for the entire 7 days. 

The laser data was assumed to be corrupted by 10 cm white noise and the data 
acquisition rate was one per minute. The only systematic errors considered 
were those due to geopotential errors. Attention was focused on those geopo- 
tential terms which are capable of producing at least a twenty -five meter per- 
turbation of the GEOS-C orbit. Zonal terms were excluded and this left twenty- 
four geopotential coefficients which were considered as error sources. The 
standard deviations of these coefficients are determined by differencing the 
values of the coefficients in the S.A.O. 69 geopotential with corresponding co- 
efficients in the G.E.M. 4 field. The significance of a geopotential term as an 
error source in state estimation is, of course, a function of the point in the or- 
bit that one considers. By utilizing mathematical procedures outlined in the 
appendix it was possible to determine at various points along the orbit the root 
sum square contribution to the uncertainty in each component of satellite state 
due to each considered geopotential coefficient. These contributions are called 
aliasing terms and the geopotential coefficients were ordered according to their 
maximum altitude aliasing along the seven day arc. The results are displayed 
in Table 1. Columns one and two contain the order and degree of the geopotential 
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Figure 1 . Visibility for GEOS-C 
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Table 1 


Maximum radial, along track, and cross track errors (in meters) of orbit de- 
termination of GEOS-C due to uncertainties of dominant geopotential terms. 


Geopotential 

Coefficient 

o 

Maximum 

Height 

a 

Maximum 

Along 

Track 

o 

Maximum 

Cross 

Track 

a 

8(6,6) 

0. 5(10)* 10 

4. 88 

19.72 

7. 16 

C(6, 5) 

0.1(10)- 9 

3.47 

14. 12 

6. 19 

C(4, 4) 

0. 2(10)' 8 

3. 0 

18. 15 

4.9 

C(6,6) 

0. 16(10)' 10 

1.88 

5. 62 

2.26 

C(3,2) 

O.l(lO)* 7 

1.4 

2.9 

0.52 

C(4,3) 

0.2(10)- 8 

1.3 

9. 0 

2.13 

8(4,3) 

0.3(10)- 8 

1.3 

11.00 

3.7 

8(5,2) 

0.3(10) - 8 

0.96 

2. 18 

0.28 

8(6,2) 

0.3(10)- 8 

0.91 

7.20 

1.20 

C(2, 2) 

0.2(10)- 7 

0. 80 

9. 9 

2.3 

8(6,5) 

0.2(10)- 10 

0. 62 

3.76 

1.23 

8(4,2) 

0. 1(10)- 7 

0. 6 

3. 6 

5.00 

8(3,1) 

0.4(10)-* 

0.43 

1.0 

0. 083 

C(3 , 3) 

0.15(10) -8 

0.37 

0. 96 

0.26 

C(5,2) 

0. 1(10)- 8 

0.29 

0.65 

0. 11 

8(3,3) 

0. 1(10) -8 

0.26 

0.64 

0.29 

C<6,2) 

0. 1(10) -8 

0.22 

2.2 

0. 38 

8(4,4) 

0. 1(10)' 9 

0.14 

0.66 

0.21 

8(3,2) 

0. l(10)-8 

0. 13 

0.27 

0. 064 

C(4, 2) 

0.3(10)-8 

0. 13 

0. 96 

1.3 

C(3, 1) 

0. l(lO)- 8 

0. 11 

0.22 

0. 026 

C(4, 1) 

0 . 1(10)- 8 

0.064 

0. 51 

0.1 

8(2 , 2) 

0. 6(10)- 9 

0. 03 

0.009 

0.08 

8(4, 1) 

0 . 1(10)- 8 

0.023 

0.23 

0. 08 
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coefficients and their associated standard deviations. Columns three, four and 
five contain the maximum aliasing contributed respectively in the radial, along 
track, and cross track direction by each term. 

It is interesting to examine the aliasing of a geopotential coefficient as it varies 
along an arc. Figures 2, 3 and 4 are plots of the contributions to the standard 
deviations of the radial, along track, and cross track components of GEOS-C 
state due to coefficient S(6, 6) for the first fifteen hours after epoch. What is 
seen at least in the radial and cross track direction is a periodic occilation 
whose period is the rectification of the orbital period and whose amplitude is 
modulated by a complex observability pattern. The pattern in the along track 
direction cannot be interpreted so simply. 

Attention is focused on the altitude determination requirements of GEOS-C im- 
posed by its on board altimeter. Since the instrument is capable of one to two 
meter precision it is required to determine GEOS-C altitude to a standard de- 
viation of one to two meters. If, say, n parameters are to be added to the esti- 
mated state in order to improve altitude determination accuracy, then the maxi- 
mum benefit from the added computational cost is obtained when the estimated 
parameters are the first n geopotential coefficients in the list of Table 1. Hence 
the cheapest way to satisfy altitude determination requirements for GEOS-C is 
to determine how many consecutive terms in the list of Table 1 must be esti- 
mated in order to meet the requirements. Using covariance techniques we de- 
termined the total standard deviation of altitude estimate for the GEOS-C space- 
craft assuming the postulated data was used to estimate state and, no geopotential 
terms, the first 3 terms of the list, the first 8 terms, the first 12 terms, and 
finally all 24 terms. The altitude errors were propagated along the entire seven 
day arc and the percentage of time that the standard deviation was under one 
meter and under two meters was computed. The results are shown in Table 2. 
The simultaneous estimation of the twelve dominant geopotential coefficients ap- 
pears to be a logical strategy for satisfying mission requirements at the small- 
est possible computation cost. A more detailed picture is provided by Figure 
5 which shows the altitude determination standard deviation for the first 15 
hours after epoch when no geopotential terms are estimated and when the domi- 
nant 12 geopotential terms are estimated. The change in the maximum standard 
deviation in the along track error is 26 meters when no geopotential terms are 
estimated and 12 meters when the dominant 12 geopotential errors are estimated. 
In the case of the maximum standard deviation in the cross track direction the 
improvement is from 10 meters to 5 meters. 

Within the context of this study it was not possible to determine the extent to 
which the estimates of geopoential coefficients obtained as a byproduct of the 
orbit determination procedure are physically meaningful. In general, geopo- 
tential coefficients of the same or of adjacent degrees are highly non orthogonal 
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Figure 2. Contribution to Standard Deviation of Radial Component of GEOS-C due to S(6,6) 
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Figure 3. Contribution to Standard Deviation of Along Track Component of GEOS-C due to S(6,6) 
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Figure 4. Contribution to Standard Deviation of Cross Track Component of GEOS-C due to S(6,6) 



Table 2 


Altitude determination accuracy for GEOS-C as a function of number of 

estimated geopotential terms. 


Number of 
Geopotential 
Terms 
Estimated 

% of time 
under 1 meter 

% of time 
under 2 meters 

0 

0 

0 

3 

8 

55 

8 

50 

100 

12 

80 

100 

24 

100 

100 


in data obtained from a single satellite. Consequently the estimated geopotential 
terms are likely to be badly aliased by uncertainties in unestimated geopotential 
coefficients. But it is precisely this aliasing effect on the estimated terms that 
permits them to be effective in functioning as lumped parameters which absorb 
the effect of uncertainties in unestimated geopotential coefficients . 
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Figure 5 . Altitude Standard Deviation for GEOS-C 



CONCLUSIONS AND RECOMMENDATIONS 


This paper has suggested that the systematic use of covariance analysis tech- 
niques is effective in identifying the cheapest estimation strategy which will 
satisfy a given set of mission requirements for orbit determination accuracy. 
The techniques were used to determine an estimation strategy which satisfied 
mission requirement of one to two meter altitude resolution for GEOS-C orbit 
determination. The procedures outlined in the paper should be useful in meet- 
ing the stringent orbit determination requirements of future applications 
satellites. 

Although the results which have been obtained are encouraging, the present 
study has several limitations. The most obvious limitation is that we consider 
only 24 geopotential coefficients as error sources for GEOS-C orbit determina- 
tion. A full field up to degree and order 8 together with certain resonant terms 
should be considered. When this is done, covariance techniques will answer the 
question of whether the estimated geopotential coefficients have physical signi- 
ficance or are functioning as lumped parameters . Another limitation is that 
covariance techniques , powerful as they are , cannot provide information about 
certain kinds of difficulties which may arise when the least squares computa- 
tional algorithm is fully implemented. The covariance results should be sup- 
ported by full simulations. The final limitation which will be mentioned relates 
to the fact that as the orbit precesses the set of geopotential coefficients whose 
uncertainties are significant error sources will change. This means that at 
certain intervals it may become necessary to repeat the entire procedure which 
is outlined in the paper. We have not considered this complication in the study. 
It is a difficulty which should be thoroughly explored in future studies. 
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APPENDIX 

COVARIANCE ANALYSIS AS APPLIED TO ORBIT DETERMINATION 


COMPUTING COVARIANCE MATRICES 

Let f(m) be an m dimensional vector consisting of the differences between the 
correct values of observations of a satellite and nominal values of the observa- 
tions as determined from a nominal orbit. Also let z(n) be an n dimensional 
vector of differences between actual and nominal values of the state of the 
satellite at an epoch and differences between actual and nominal values of 
parameters in the dynamic and measurement models whose associated un- 
certainties may limit our ability to estimate satellite state from the data. The 
sensitivity matrix c (m, n) is defined as that matrix whose element in the ith 
row and the jth column is the partial derivative of y(i) with respect to z(j). A 
first order Taylor series expansion of the functional relationship between y and 
z about the nominal value of z yields 


y 


'V. 

C Z 


(A-l) 


An orbit determination program in processing observations y of y to obtain a 
least square adjustment to z computes a so-called normal matrix defined as 

r , (n, n) - c T w c (A-2) 


where w is a weighting matrix and is usually the inverse of the covariance 
matrix of the observations y of Once an orbit determination program com- 
putes and stores the normal matrix, a number of questions can be raised and 
answered at very little cost in terms of computation time. 

The best estimate of the state of the satellite at epoch is obtained by perform- 
ing a least squares adjustment of the state at epoch and all other parameters 
with which are associated significant uncertainties. But frequently this straight- 
forward approach leads to severe core storage requirements. In practice some 
of the parameters in the dynamic and measurement models are estimated along 
with state and others are fixed at their nominal values and left unadjusted in the 
least squares process. In order to determine the consequences of estimating 
some parameters and ignoring others it is useful to compute the covariance 
matrix of such a least squares estimation procedure. 
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Let 'i be decomposed into two disjoint parameter sets as follows 


'V* 

z 


x,(nj) 


^ 2 (n 2 ) 


(A- 3 ) 


where x p is a set of n lf parameters which are to be estimated in a least squares 
process and x 2 is a set of n 2 parameters whose nominal values are left un- 
adjusted by the least squares process but whose uncertainties are to be con- 
sidered in computing the covariance matrix of the resulting estimator. Define 
a matrix A(m, q,) as a matrix whose element in the ith row and jth column is 
the partial derivative of y(i) with respect to x t (j). Analogously define B(m, n 2 ) 
as the matrix whose element in the ith row' and jth column is the partial deriva- 
tive of y(i) with respect to x 2 (j). For future reference notice that the normal 
matrix 17 of z as computed and stored by an orbit determination program and 
defined by Equation A-2 can be written as 


V 


~ A T wA A T wB " 
_B T wA B t wB_ 


(A- 4 ) 


Assume that there exists a priori estimates of Xj and x 2 with properties 


/ _ ~ 

Xj - Xj + 04 , 

E(a x ) = 

Ol 

E(a a T ) 

1 

X 2 X 2 CE. 2 » 

E(a 2 ) - 

0, 

E ^ cl 2 a 2 ) 


and assume that the observation vector y or y has properties 


y - y + v, E(v) - 0 , E(vi' T ) = w 1 


The least squares estimate of Xj is obtained as the value of xi which minimizes 
the loss function 


L ( x i) - (y - Axi - Bx 2 ) T W(y - Axj - Bx 2 ) + (Xj - Xj ) 1 Pj 1 (x/ - Xj) (A-5) 
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The resulting least squares estimator of x 1 is well known to be 

x x = (A T wA + P" 1 )” 1 [A T w(y - Bx a ') + Pf 1 x t '] (A-6) 

Define 

P = [E (x x - x x ) (x, - x x ) T ] (A- 7) 

A series of substitutions reveals that 

Xj - Xj = (A t wA + P^ 1 ) -1 (- A T wBa 2 + A T wv + Pj 1 a.j) (A-8) 


Equation 8 yields 

P = (A t wA + P" 1 )' 1 + (A T wA + P; 1 )" 1 A t wBP 2 B t wA(A t wA + Pj -1 )" 1 . (A-9) 

Notice that the right side of Equation 9 can be computed if one has a priori co- 
variance matrices P t and P 2 , and the upper right and upper left portions of the 
normal matrix. To determine the covariance matrix of an estimator which 
estimates some subset of 2 other than x : , all that is necessary is to permute 
the rows and columns of 17 in the appropriate fashion and proceed as before. 
Thus if one assumes that the normal matrix defined by Equation 2 is precom- 
puted it becomes an easy matter to obtain the resultant covariance matrix when 
any subset of the z parameters are estimated in a least squares sense and the 
rest are ignored. 


THE ALIAS MATRIX 

Assume that all the data has the same variance. Hence 

w - (Icr 0 2 >- 1 (A- 10) 


where cr 0 2 is the common variance of each data point. Also assume that the a 
priori estimates of the unadjusted parameters are independent. Under this 
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as sumption the covariance matrix P, of x L can be written as 


P 2 ^ 


0 


(A- 11) 


where o-. 2 is the a priori variance of the ith unadjusted parameter. Also define 
a matrix K(ni, n 2 ) as 


K, = ( A t wA) - 1 A T wB 


(A- 12) 


With these assumptions Equation 9 yields the following expression for the ith 
diagonal element of P 


n 2 

p(i ’ i) = (/?i 'j CT j )2 (A-i3) 

i=o 

where /2. 0 <is the ith diagonal element of the matrix (A T A) -1 (this assumes that 
diagonal elements of the matrix P^ 1 are relatively small) and 


0 itj - K(i,j), j>l (A-14) 

The standard deviation of the ith estimated parameter is given by 


cr. 


1 



Define the error sensitivity matrix as 


(A -15) 


S = 




= 1, 2, 


j = 0, 1, 


n, (A- 16) 
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And finally define the Alias Matrix as 


L = Sa 


(A- 17) 


where 


cr 



0 


0 



(A- 18) 


The standard deviation of the ith estimated parameter is seen to be the root sum 
square of the terms in the ith row of the alias matrix. The elements in the first 
column of the alias matrix represent the RSS contribution to the standard devia- 
tion of each estimated parameter due to the data noise. The elements in the jth 
column, j > 2, represent the RSS contribution to the standard deviation of each 
estimated parameter due to the j - 1st unadjusted parameter. 

Possession of the alias matrix reveals much of the probability structure of the 
postulated least squares estimator. With this information one can quickly de- 
termine which error sources are significant with regard to the estimation of a 
given parameter. 

Propagating Covariance Matrices 


Equation 9 provides the covariance matrix of the state x 1 at some specified 
epoch. In many cases it is important to determine how accurately the state 
can be determined at some time other than epoch. In order to do this cor- 
rectly it is necessary to take into proper account uncertainties in dynamic 
parameters. These parameters may be in an estimated mode or in an un- 
adjusted mode and to incorporate their effect one resorts to state transition 
matrices which presumably have been precomputed by an orbit determination 
program. Let x x (T) be the estimated state at time T. Assume as output from 
an orbit determination program the state transition matrices 


v x (T) 


3*i (T) 



v 2 (T) 


dx, (T) 
<?x 2 


(A- 19) 
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If there are no dynamic parameters in the estimation vector x lt the matrix 
Vj (T) takes on the particularly simple form, 


v x (T) 


S 0 

0 I 


(A-20) 


where S is the six by six matrix defined as the partial derivative matrix of the 
state of the satellite at time T with respect to the state of the satellite at epoch. 
If dynamic parameters are included in the estimated state, the off diagonal 
matrices become non- zero and Vj (T) assumes a more complicated form. 

The matrix v 2 (T) is the matrix of partial derivatives of the state (T) with 
respect to the unadjusted parameters x 2 . If no dynamic parameters are in 
the unadjusted mode, v 2 (T) is the null matrix. A first order Taylor series 
expansion of the function which describes the time evolution of the state x (T) 
yields 1 


*i ( T ) " v i ( T ) *i + v 2 (T) x 2 (A- 21) 

Substituting as obtained from Equation 6 for x 1 and x 2 ' for x 2 provides the 
best estimate x t (T), of x 1 (T) 

x i ( T ) = v, (T) xj + v 2 (T) x 2 ' (A-22) 

The covariance matrix of x 1 (T) is given by 

p (T) =. v, (T) PvJ (T) + v 2 (T) P 2 v 2 t (T) + v, (T) E[x 1 x‘] v 2 (T) 

+ v 2 (T) E [x ' xjl vj (T) (A _ 23) 

Equation 23 in conjunction with Equations 6 and 9 yields 

P (T) = v t (T) (A T w A + P " 1 )" 1 Vj T (T) +[ v, (T) (A T w A + P " 1 )" 1 A T wB 

- v 2 (T)] P a [v x (T) (A t wA + P - 1 )" 1 A T wB - v 2 (T)] T (A-24) 
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Finally notice that in much the same fashion that Equation 9 was used to de- 
velop an alias matrix at epoch, Equation 24 can be utilized to develop an alias 
matrix for any time T. 


REMARKS 

If one possesses a functioning orbit determination program it becomes a rela- 
tively easy matter to add covariance analysis capability to the system, A com- 
puter program can be written which assumes as input a normal matrix and state 
transition matrices as generated by the orbit determination program. By 
permuting the rows and columns of the normal matrix and completing the 
matrix operations defined by Equation 9, the covariance matrix of a least 
square process which adjusts any subset of the parameters and ignores the 
rest can be computed. An alias matrix can be obtained and significant error 
sources can be identified. By utilizing the precomputed state transition matri- 
cies, the covariance matrix of the estimate of the state can be propagated from 
epoch to any other time. These operations are very simple and they consume 
little computer time. 

Since the normal matrix and state transition matrices are computed once and 
permanently stored, it is possible to investigate a large number of possible 
estimation strategies. This can be done conveniently and cheaply. For many 
applications such a program is a useful and quickly developed addition to an orbit 
determination system. 
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